Resonant peak in the density of states in the normal metal / diffusive ferromagnet / 

superconductor junctions 
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The conditions for the formation of zero-energy peak in the density of states (DOS) in the nor- 
mal metal / insulator / diffusive ferromagnet / insulator / s-wave superconductor (N/I/DF/I/S) 
junctions are studied by solving the Usadel equations. The DOS of the DF is calculated in various 
regimes for different magnitudes of the resistance, Thouless energy and the exchange field of the 
DF, as well as for various resistances of the insulating barriers. The conditions for the DOS peak 
are formulated for the cases of weak proximity effect (large resistance of the DF/S interface) and 
strong proximity effect (small resistance of the DF/S interface). 



In ferromagnet/superconductor (F/S) junctions 
Cooper pairs penetrating into the F layer from the S 
layer have a nonzero momentum due to the influence of 
exchange fieldiSi^. This results in oscillating behavior 
of the pair amplitude or a 7r-phase shift of the order 
parameter in the ferromagnet. A negative sign of the 
real part of the order parameter may occur when the 
thickness of the F layer is larger than the coherence 
length of the F layer. The occurrence of the 7r-phase shift 
makes it possible to realize the SFS tt -junctionsi, as was 
confirmed experimentally^*^*^*^. The order parameter 
oscillations also lead to nonmonotonous dependence of 
Tc in SF bilayers on the F-layer thicknes92ii2iiiii2J^. 
Effects of resonant transmission in conductivity of SF 
structures were discussed in Refl^*i^*i&. 

Another interesting consequence of the oscillations of 
the pair amplitude is the spatially damped oscillating 
behavior of the density of states (DOS) in a ferromag- 
net predicted theoreticallyiliiSiiSiSli in various regimes. 
The energy dependent DOS calculated in the clearilS 
and the dirtySi limits exhibits rich structures. Exper- 
imentally DOS in F/S bilayers was measured by Kon- 
tos et al. who found a broad DOS peak around zero 
energy when the 7r-phase shift occursSS. In diffusive 
ferromagnet/superconductor (DF/S) junctions the zero- 
energy DOS may have a sharp peak^^. However the con- 
ditions for the appearance of such anomaly have not been 
studied systematically so far. 

The purpose of the present paper is to calculate DOS 
in N/DF/S junctions and to formulate the conditions for 
the zero-energy DOS peak in two regimes correspond- 
ing to the weak proximity effect (large DF/S interface 
resistance) and strong proximity effect (small DF/S in- 
terface resistance). We will show that in the former case 
the condition is equivalent to the one of Refill, while in 
the latter case the new condition is found. The calcu- 
lation will be performed in the zero-temperature regime 
by varying the interface resistances as well as the resis- 
tance, the exchange field and the Thouless energy of the 
DF layer. 

We consider a junction consisting of normal and 
superconducting reservoirs connected by a quasi-one- 



dimensional diffusive ferromagnet conductor (DF) with a 
resistance Rd and a length L much larger than the mean 
free path. The DF/N interface located at x = has the 
resistance i?^, while the DF/S interface located at a; = L 
has the resistance Rb- We model infinitely narrow insu- 
lating barriers by the delta function U{x) = HS{x — L)-\- 
H'S{x). The resulting transparencies of the junctions T„i 
and are given by Tm = 4cos^ (/)/(4cos^ (j) -\- Z'^) and 
T/„ = 4cos2(/)/(4cos2 Z'^), where Z = 2H/vf and 
Z' = 2H'/vp are dimensionless constants and cj) is the 
injection angle measured from the interface normal to 
the junction and vp is Fermi velocity. 

In the following calculation we will apply the quasi- 
classical Green's functions formalism. The 2 x 2 re- 
tarded Green's functions in N, DF and S are denoted by 
Rq{x), Ri{x) and R2{x) respectively. Ro{x) and R2{x) 
are expressed by Ro{x) = and R2{x) — {gr^ + /f2) re- 
spectively, with g = e/y/e^ — and / — A/VA^ — e^, 
where T2 and fa are the Pauli matrices, and A and s de- 
note the energy gap and the quasiparticle energy mea- 
sured from the Fermi energy, respectively. It is con- 
venient to use the standard -parametrization when 
function Ri{x) is expressed as Ri{x) = t^cosO^x) + 
T2 smOix). The parameter 9{x) is a measure of the prox- 
imity effect in DF. The spatial dependence of 9{x) in DF 
is determined by the static Usadel equationS^ 



D 



—I 

dx"^ 



{x) + 2i{e - {+)h) sin[6'(a;)] = 



(1) 



for majority (minority) spin with the diffusion constant 
D and the exchange field h in DF. Note that we assume a 
weak ferromagnet and neglect the difference of the Fermi 
velocities of the majority and minority spin subbands. 

Further we shall apply the Nazarov's boundary 
conditionals^ for e{x) at both interfaces. At the DF/N 
interface it has the following form: 

<F >' 

(2) 



L 80 (x) 
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(2-T;J+r^cos0(O+)- 
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and it has a similar form at the DF/S interface. This 
boundary condition is based on the Zaitsev's bound- 
ary conditions^ with isotropic hmit and generahzes the 
Kupriyanov-Lukichev boundary conditionSl. 

The average over the various angles of injected parti- 
cles at the interface is defined as 



< Bi^) >'= 



cos 



d^T '{4,) COS cl, 



with B{(l)) = B and T'((/)) = T/^. The resistance of the 
interface i?^'-* is given by 



n(') _ r>i') 
- ^0 



j:i%dcj,T('){ct,) COS cj,' 



Here, for example, i?^''' denotes i?f, or i?^, and Rq' is 
Sharvin resistance, which in three-dimensional case is 
given by i?Q ' ^ = e^kpSi'^ /{An'^), where kp is the Fermi 

wave-vector and Sc'' is the constriction area. 

In the following, we will study the local DOS N in the 
DF layer which is given by 



?(') 



N/No = i^Recos6'(a;) 



hi 

where A^o denotes the DOS in the normal state. The DOS 
will be calculated by numerical solution of the Usadel 
equations with the boundary conditions given above. 

Below we will concentrate on the DOS at a; = (N/DF 
interface) in the regime of large resistance of the N/DF 
interface, Rd/R'j^ <C 1 and will also fix the barrier trans- 
parency parameters Z = 3, Z' = 3. 

In order to study the condition for the appearance of 
the zero energy DOS peak, we plot the normalized zero 
energy DOS at x = as a function of Eth ~ Djl?. 
Fig. n shows the DOS for Rd/R'b = 0.1 and various 
/i/A. In Fig. n (a) the zero-energy peak appears at 
Exh ^ "^hRb/Rd, while in Fig. H (b) and (c) the peak 
appears at Exh ~ h. Thus we can conclude that the 
condition for the DOS peak for large Rd/Rb is essentially 
different from the one for small Rd/Rb- 

Fig. 13 shows the DOS as a function of e for the pa- 
rameters corresponding to the peaks in Fig. ^for various 
h/A. In all these cases the DOS peak appears around 
zero energy. For small h/A the DOS peak is narrow but 
it becomes broader with the increase of h/A. It's im- 
portant to note that this peak not always requires the 
sign change of pair amplitude. This is also clear from 
the fact that the peak occurs for large Thouless energy 
(short DF) when there is no sign change. For other set of 
parameters the DOS peak is smeared as they break the 
condition Exh ~ '^hRb/Rd or Exh ~ h. 

Let us first discuss the case of strong proximity effect 
in detail. Fig. |21 shows the zero energy DOS at x = as 
a function of Exh for h/A = 1 and various Rd/R'i, with 
(a) Rd/Rb — 5 and (b) Rd/Rb — 10. In this case the 
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FIG. 1: Normalized zero energy DOS as a function of Exh 
for large resistance of the N/DF interface Rd/R'i, = 0.1 and 
various h/A with resistance ratios at the DF/S interface (a) 
Rd/Rt = 1, (b) Rd/Rb = 5 and (c) Ra/Rt = 10. 



peak at Exh ^ h is suppressed with increasing Rd/R'f,. 
Therefore this condition is valid for small Rd/R'i^ 

Fig. 21 shows the spatial dependence of Im6' for ma- 
jority spin for Rd/R[, = 0.1, Exh/A — 1 and vari- 
ous h/A with (a) Rd/Rb = 5 and (b) Rd/Rb = 10. 
For the appearance of the DOS peak, large value of 
Ini0 is needed because the normalized DOS is given by 
Recos(6') = cos(i?e(6')) cosh{Im{9)). As seen from Fig. 01 
the magnitude of lm6 increases with the increase of the 
distance from the DF/S interface and achieves a maxi- 
mum when Exh — h. 

Note that the zero-energy DOS at a; = does not de- 
pend on Exh if the condition Exh = h holds. To explain 
that, let's write Eqs. ^andl^lat £ = 0: 
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FIG. 2: Normalized DOS as a function of e for i?d/7?6 = 0.1 
and various h/A with (a) Rd/Rb = 1 and _ETh = 2hRb/Rd = 
2/1, (b) Rd/Rb = 5 and _ETh ~ h, and (c) Rd/Rb ~ W and 
^Th = h. 



where y = x/y/D/h. Since for Eth = h we have 
D/h = ExhL'^/h = L^, the above equations don't con- 
tain Eth as a parameter. Similar arguments can be ap- 
plied to another boundary condition at DF/S interface. 
This proves the above statement about independence of 
the zero-energy DOS at a; = on Eth- 
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FIG. 3: Normalized DOS at zero energy as a function of Eth 
for /i/A = 1 and various Rd/R'b with (a) Rd/Rb = 5 and (b) 
Rd/Rb = 10. 



form: 

cos6'o| = 



Ri 



2i[e-h) 
Eth 



(3) 



Rd I Rd„ 
Ri ^ Ri^ 



2i{e-h) 



For Rd/R'ij = and e = 0, the DOS has the form 



cos Uqi 



2ih 



( Ea\ - ( JUL. 

\Rb J \ Eth 



(4) 



Now let us discuss the weak proximity effect and derive 
the condition Rd/Rb "^^h/ETh, following Ref.'^. When 
spatial variation of 6 is small, i.e., L <C \/E>/ \ h \ (for 
the spin-up or spin-down subband respectively) and both 

Rd/Rb and Rd/R'b are small (weak proximity effect), 9 Eg ^ ExhRd/Rb for weak proximity effect, or Eg ^ E' 
can be expanded as6 — 6o + 9ix + 92X^ where 0i, 6*2 ^ ^o- 
Note that the derivatives of 9 are proportional to these 
quantities at the interfaces (see Eq. (0) and RefiSi). 



which provides the resonant condition Rd/Rb ^ '^h/Exh- 
Similar result follows for the spin-down subband by re- 
placing h by —h. 

Another resonant condition for the strong proximity ef- 
fect, Eth ~ h, is equivalent to the condition L ^ \J D/h. 
Thus, zero-energy DOS peak appears when the proximity 
effect is strong and the length of ferromagnet is of the or- 
der of the coherence length in a ferromagnet S^p — yD/h. 

Let us discuss the physical meaning of two condi- 
tions. In DN/S junctions there is a minigap Eg^ where 



Th 



for strong proximity effect^^. In DF/S junctions this 
minigap is shifted by h, then the DOS peak appears when 



In this case the solution of the Usadel equation in the 
spin-up subband satisfying boundary conditions has the 



h^Eg. 

Note that in the calculations we have fixed Z = Z' = ?>, 
but the specific parameter choice does not change the 
results qualitatively. 
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In summary, we have studied the conditions for the 
appearance of the DOS peak in diffusive ferromagnet, 
in normal metal / diffusive ferromagnet / s-wave super- 
conductor junctions. We have discussed two regimes of 
weak and strong proximity effect depending on the ra- 
tio Rd/Rb- The results in the regime of weak proximity 
effect are essentially the same as found in Refi^. How- 
ever, in the regime of strong proximity effect the results 
are qualitatively different. Let us summarize the two 
conditions: 

1. When the proximity effect is weak [Rd/Rb ^ f), 
the condition for the DOS peak is Rd/Rb ^ 2h/ETh- 

2. When the proximity effect is strong (Rd/Rb ^ 1), 
the DOS peak appears when Eth ~ h, i.e. when the 
length of ferromagnet is of the order of the coherence 
length y/D/h. 

Note that the above two conditions cross over into each 
other when Rd/Rb ^ 2. Since the DOS is a fundamental 
quantity affecting various physical properties, our results 
may have many applications, e.g., for the conductance of 
N/DF/S structures. 



FIG. 4: Spatial dependence of lm6 for majority spin for 
Rd/R'b = 0.1, Eth/ A = 1 and various h/A with (a) Rd/Rb = 
5 and (b) Rd/Rb = 10. The DF/N interface and the DF/S 
interface are located at a; = and x = L respectively. 
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